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IMPROVED DEFINITION OF CRUSTAL ANOMALIES FOR MAOSAT DATA 


ABSTRACT 

MAOSAT vactor magnatomatar data can ha routinalg corractad for axtarnal 
flald affacta auch aa tha ring currant and tha dally variatlonf by 
flltarlng long-wavalangth harmonica from tha data. Saparation of 
flalda dua to low**alt Ituda aourcaa from thoaa of hlgh**altltuda aourcaa 
1* affactad by maana of dual harmonic aipanalona In tha aolution of 
Dlrichlat'a problam. Thla ragraaalon/harmonic filtar procadura la 
appliad In an orbit-by-orbit baaia« and Initial taata on MAOSAT data 
from orbit 1176 ahow raductlon in axtarnal flald raalduala by 24.33 nT 
RMS in tha horiiontal componant* and 10.93 nT RMS in tha radial 
componant. 

INTRODUCTION 

Tha af f ac t i vanaai of tha MAOSAT miaaion for dallnaatlon of cruatal 
magnatic anomaliaa ia hamparad by tha praaanca of axtarnal fiald 
contamination in tha data. Flalda from tha ring currant and tha dally 
variation ara aa larga aa that of all but tha largaat cruatal magnatic 
anomaliaa at aatalllta altituda* (Sugulra and Hagan* 1979). Tha 
pattarna of tha ring currant and tha dally variation ara hamiapharlcal 
or global in acala* but bacauaa thay can changa auddanly In tlma* tha 
aaaoclatad axtarnal flalda maaaurad at tha aatalllta can aliaa with 
thoaa of tha ganarally ahortar uavalangth cruatal anomaliaa. Whila tha 
groaa faaturaa of tha ring currant and tha dally variation ara fairly 
wall known* thay ara not amanabla to pradlctlva phyalcal modallng 
bacauaa of thair unpradictabla variability in raaponaa to tha aolar 
wind. Suglura and Hagan C19791 hava found that tha Solar guiat (Sq) 
daily variation It wall rapraaantad by apharical harmonica of dagraa 4* 
for axampla* but thaaa modal coafficlanta vary from hour to hour. Far 
from balng a pradlctlva modal* thla rapraaantatlon ia puraly 
daacriptlva* balng baaad on magnatograma from many magnatic 
obaarvator iaa dlatributad around tha globa. An axampla of tha affacta 
of thaaa lonoapharic curranta on aatalllta data ia illuatratad in 
Figura 1* ahowing aavaral PODO aatalllta paaaaa in which tha fiald dua 
to tha dally variation ia apparant and la varlflabla by comparlaon with 
obaarvatory data. Although thla data la far ramovad in local tlma from 
tha noon-tima paak in daily-variation amplltuda. Ita affact la atlll 
aignificant comparad to that of tha cruatal magnatic anomaliaa of 
intaraat. Plata of tha flald dua to tha Sq. currant ayatam on Dacambar 
27* 1964* 1700 hra OMT ara ahown In Figura 2 for varleua latltudaa and 
altltudaa. Not^a that thla tlma corraaponda to a parlod of minimum 
aolar activity*- whila tha MAOSAT miaaion waa flown during a aolar 
maximum. Figura 3 illuatrataa tha incraaaa In Sii affacta to ba 
axpactad in a parlod of high aolar activity. 

• ring currant ia aomawhat aaaiar to modal* alnca tha gaomatry ia 
almplar and tha Oat indax daacribaa Ita atrangth. Figura 4 ahowa fiald 
raaiduala for aavaral paaaaa ovar tha aama araa for diffarant aolar 
wind condltlona* aa avldancad by tha varlabllllty of tha Oat indax. 
Tha aolid linaa in aach plot ara tha corraction daviaad by Cain and 


Paga 


2 


PHOENIX CORPORATION 


March 19> 1982 


Davis tl9733 for this offset. This corroction is basod on a loast 
s«iuaros curvo fitting to the obsswod data# a procodura which raiiuirss 
caraful attantion to tha apoch and longth of tha data arc to minimiia 
aliasing with tha crustal signals of intarast. A physical modal of tha 
ring currant* with currant dansity varying as a function of tha Dst 
indicas* is attractiva bacauss of its simplicity. Such a modal ramains 
basically dascriptiva in naturo.^ howavor* baing dapondont on tha 
analysis of obsarvatory magnatograms for tha Dst indicas* which may 
changa hourly. 

Othar procaduras for ramoving axtarnal fiald offsets from tha 
satallita magnatomatar data ara basad on minimizing tha variability 
batwaan individual data profilat ovar tha sama araa t difforont timas. 
Raagan* at al. C19793* usod a sami-'quantitativa visual intorcomparison 
tachniqua for aditing savaraly contaminatad passas. Mayhaw C19793* has 
omployad tha fitting of linaar or quadratic functions of varying 
longths to minimita tha intarpass variability. 

Tha problom with tha above techniiiuas is that thay ara all somawhat 
ad hoc in natura* and ara not easily gonarallzed into an algorithm for 
routine processing of satellite magnatomatar data. In general* the 
data must bo carefully salactad by geographic area* arc length* and 
time* and a year or more may elapse before enough observatory data can 
be collected and analyzed to produce models of 8g or Dst indicas. Evan 
at bast* tha resulting models of tha daily variation and tha ring 
currant ara global averages* and cannot describe tha sudden variations 
in time and space which can occur in these currant systems. 
Furthermore* these models have been developed to deal only with scalar 
total fiald data* not the more variable vector fiald components 
provided by tha HACSAT mission. 

This report describes the results of an investigation to develop a 
method for routine removal of tha most persistant axtarnal fiald 
affects from tha MAOSAT data. Tha proposau approach involved modifying 
existing ring currant correction procaduras to use indices such as DS 
and Dst* and simplifying tha Sugiura and Hagan modal for Sg currants. 
Lack of timely data from magnetic observatories* and early success with 
raoression models for correcting the ring current* led us to adopt a 
different approach. A prototype method which uses knowledge of tha 
spectral content of the external fields* but which is insensitive to 
geographic araa and can be used routinely in 
naar-to-tha'-maasuramant'time* was developed. Tha following sections 
describe this prototype method and tha results of its application to 
MADSAT vector magnetometer data. 

EXTeRNAL FIEUP ..wPEL 

In tha early phases of this investigation* an evaluation was made to 
assess tha application of existing models for tha daily-variation and 
tha ring currant. These efforts concentrated on tha addition of Dst 
indicas to Cain's ring currant modal* and on simplifying tha Sugiura 
and Hagan C19673 modal for tha daily variation to a procedure more 
amenable to prompt* routine processing. 
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aina.fiyrrtnt 

To fully appraciata tha maaning of tha Dat indax. It it iaportant to 
undaratand Ita darivatlon. Tha ring currant la thought of aa a aiora or 
laaa ayaaiatric ring cf currant at about 5 aarth radii, locatad in tha 
gaoaagnatic aquatorial plana. This currant ia ganaratad by tha drift 
of protona and alactrona trappad baloui 7.9 aarth radii batwaan tha 
dipolar main fiald and tha magnatopauaa. Tha drift ia a function of 
tha radial gradiant in magnatic fiald intansity. and raaulta in a nat 

waatward currant vactor. During magnatic atorma. tha incraaaad 

injaction of particlaa into tha magnatoaphara. and tha incraaaad aolar 
wind praaaura cauaa tha ring currant to intanaify and bacoma loan 
aymmatric. dapraaaing tha magnatic fiald atrangth obaarvad at tha 

oarth'a aurfaca by aa much aa 100 nT. Sugiura C19643 da/alopad tha Dat 
indax aa an indicator of thi. «trangth of tha ring currant. Ha modola 
tha fiald dua to tha ring currant aa a Fouriar aoriaa in tha 

gaomagnatic aquatorial plana aat 

D = (T)sin [(T + X) + aj^(T)] + Aq(T) (1) 


whara T ia univaraal tima. X ia local tima. Aq. Ax« nnd a 
ara modal coaff icianta. corraaponding to amplitudaa and phaaa 
anglaa. raapact ivaly. 

By obaarving tha horizontal fiald componant at 4 tarraatrial 
obaarvator iaa. and adjuating tha obaarvationa for tha latituda of tha 
atation. ona may aolva for thaaa amplitudaa and phaaaa. Tha Oat indax 
ia A(). tha aymmatrical fiald atrangth dua to tha undiaturbad ring 
currant at tha gaomagnatic aquator. Tha amplituda Ax ia tha D8 indax. 
rapraaanting tha aaymmatric componant of tha ring currant during 
diaturbancaa. Excapt for thaaa aaaymatriaa in tha ring currant fiald. 
it ia approx imataly uniform and parallal to tha aarth 'a dipola axia in 
tha vicinity of tha aarth. Thua it aliaaaa with othar fialds of 
aimilar charactar iatica. auch aa thoaa of tha magnatotail and 
magnatopauaa currants. 

A pralifflinary avaluation waa parformad of tha Cain and Davia C19733 
ring currant corraction in compariaon with tha uaa of boat fitting 2nd 
dagraa polynomials, bast-fitting 3rd dagraa cosina aariaa. and a modal 
using tha first thraa tarms of tha apharical harmonic axpanaion for a 
uniform fiald with coafficiants proportional to tha Dat and D8 indicaa. 
6avan colinaar FOOD data paaaaa wara aalactad having continuous 
covaraga from 50 dagraaa South to 50 dagraaa North ovar aach of two 
taat araas. Thasa data wara corractad using aach of tha four 
tachniquas daacribad abova. All of tha tachniquaa gava vary aimilar 
appaaring raaulta. with no significant diffarancoa in RMS of raaiduala 
or commonality of paasas. Tha total fiald dua to tha ring currant 
corraction waa roughly proportional to coa2^. whara a ia gaomagnatic 
latituda. 
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Solar Quiat Dailu Variation 


In dava loping a modal for Sq* Sugiura and Hagan C 19673 uaad a 
apharical harmonic formulation for tha aaaociatad acalar potantial 
which accountad for tha dlract 8q fiald and Ita aacondarg inducad fiald 
in tha aarth for a tarrratrial obaarvatorg. Tha acalar potantial ia 
writtan aa: “ n / r «(«) «(e) 

(T) coanX 4- h (T) aln 
n 


W(T,r,e,X) - a 


* n 
g 

_ n 


n«l 

«(i) 

(T) cosaX 


[•: 


• F V VII m 


■( 1 ) 
h (T) sinaX 
n 


] r 


P (cv^) 
n 


( 2 ) 


whara r ia radial diatanca of tha atation from tha gaocantar* 
6 ia gaographic colatituda* X tha gaographic longituda* a ia 
tha maan radiua of tha aarthi P* <coa3> ia tha aaaociatad 
Lagandra function of dagraa n and ordar m< and g” < hj • 
auparacriptad bg a or i for axtarnal and intarnal” ara tha 
coafficianta to ba datarninad aa functiona of tha univaraal 
tima T. 


Tha acalar magnatic fiald intanaitg< darivad from thia axpraaaion bg 
application of tha gradiant oparator< ia maaaurad hourlg at about 60 
ground atationa diatributad around tha world« for a pariod of daga. 
Thaaa obaarvationa ara uaad in a laaat aq,uaraa aolution for tha 
coafficianta g-j^ and h^ up to dagraa and ordar 4. A naw aolution ia 
mada for aach houra data. Whila tha 4th ordar modal coafficianta ara 
found to ba adaguata to rapraaant tha dailg variation for a givan 
houri Sugiura and Hagan C19793 find conaidarabla tranaiant variabilitg 
in tha currant pattarn from hour to hour. Indaad« Davidaon and 
Hairtslar C19683 auggaat that tha rapid variation background in 
tarraatrial obaarvationa of tha magnatic fiald ia dua to turbulanca in 
tha 8g currant agatam. 

In our avaluation of thia modal* wa initiallg concantratad on 
almplifging tha 48 tarm apharical harmonic modal to ona with fawar 
tarma and with coafficianta proportional to raadilg availabla indlcaa 
of tha atata of tha aolar wind* auch aa Dat* Kp* or Ap. Raaulta for 
P080 aatallita data wara Inconclualva* dua to tha praaanca of orbit 
arrora. Our afforta than turnad to coilacting* procaaaing* and 
analgiing magnatograma from magnatic obaarvatoriaa for tha MAOSAT 
miaaion tlma pariod In ordar to conatruct a rafaranca modal for tha 8q 
curranta. 


EYtlyitlBD. SMiMury 

Tha apatial variation of tha ring currant ia aufficlantlg aimpla 
that a ragraaaion analgala on tha actual aatallita data uaing a low 
ordar harmonic modal la found to parform aa wall aa mora aophiatlcatad 
modala uaing praacribad 08 and Dat indicaa. 8inca thaaa indicaa muat 
ba davalopad from analgaia of magnatic obaarvatorg racorda* and ainca 
thag rapraaant global avaragaa of tha ring currant ovar hour-long tima 
intarvala* tha ragraaaion analgaia approach ia mora auitabla for 
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routina corraction of tatallita magnatomatar data. 

Thaaa aama argumanta apply to cha routina corraction of tha daily 
variation. Aftar naa^ly two yaari of collacting and procaaiing 
aiagnatic obaarvatory magnatograoiai wa atill hava inaufficiant data to 
Modal tha ionoapharic 8q curranta for tha HA08AT tiaia fraaia. Efforta 
to raduca tha numbar of tarma in tha praaantly availabla Sugiura and 
Hagan modal, and to davalop coafficianta which wara proportional to 
aagnatic indicaa yialdad raaulta which wara inconcluaiva in taating 
with POOO data. Thia approach atill doaa not provlda for corractiona 
in caaa of auddan changaa in tima or apaca of tha 8q currant 
diatr ibution. 

Our aaarch for a aimplifiad mathod which can ba uaad in a routina 
mannar turr«d to uaing tha Sugiura and Hagan mathamatical formulation 
in a ragraaaion moda. Ua alraady know tha harmonic contant of tha 
daily variation fiald wall anough to accurataly modal its spactral 
bahavior in satallita data. Such an approach would ra<iuira only tha 
aiistanca of tha satallita data and ba battar abla to copa with auddan 
changas in tha 8q currant distribution. 

tfPBQACH 

It is apparant from tha foragoing raviaw that both tha ring currant 
and tha daily variation can ba raprasantad by a rathar low ordar. long 
wavalangth. modal, albait ona which might changa rapidly in spaca and 
tima. But rathar than craata a saparata corraction modal for aach 
phanomana. wa can considar a singla harmonic modal to corract for all 
long wavalangth axtarnal fialds manifast in an orbital arc of HAOSAT 
data by ragrassion analysis on tha data. Tha mathamatical statamant of 
this approach involvas racasting aquation 2 in tha form of Dirichlat's 
problam at tha satallita orbital altltuda. 

If wa know tha valua of a harmonic function on a boundary surfaca. 
than wa may computa tha valua of tha function avarywhara it ramains 
analytic. In tha casa of aquation 2. tha boundary surfaca is tha 
surfaca of tha aarth. Tha axtarnal functions, baing analytic abova tha 
aarth. dascrlba tha incidant axtarnal fiald and its raflaction from tha 
aarth's surfaca. Tha intarnal functions, analytic insida tha aarth. 
dascrlba tha fialds from tha currants inducad by tha daily variation. 
In our application. tha boundary it tha MAGSAT orbit, which wa assuma 
is a closad curva in inartial spaca. and tha axtarnal functions and 
intarnal functions maraly dascrlba tha long wavalangth bahavior of tha 
scalar magnatic potantial abova and balow tha orbit altituda. 
raspactivaly. Thus thasa functions, which ua shall tarm innar and 
outar functions to avoid confusion with tha convantional dafinition of 
tha intarnal and axtarnal fiald. raprasant tha fialds from all sourcas. 
ragardlass of thair sourca location. Tha 8q daily variation, tha 
fialds of tha ring currant, inducad currants in tha aarth. and tha 
convantional intarnal fialds of tha crustal anomalias and tha main 
dipola. ara raprasantad abova and balow tha orbit by tha outar and 
innar functions, raspactivaly. Wa affactlvaly ramova tha contribution 
of tha main dipola and tha crustal anomalias to tha long wavalangth 
fialds in thasa ragions by subtracting tha scalar potantial of tha main 
dipola and a rafaranca fiald such as tha HOST 6/80 (Langal at al. . 
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1981). With this subtraction, and ths substitution of aaan orbital 
radius for tho aoan oarth radius, aquation 8 doscribos tho rosidual 
scalar potontial for tho long wavolangth oitsrnal fialds ouch as tho 
ring currant and tba daily variation, 

8inca tha aiodal rograssion fit is to ba appliad to a singla orbit of 
data, tha spharical haraionic aipansion of aquation 2 mag ba roplacod bg 
a olanar harmonic aipansion of tha scalar magnatic potantial as: 

N - n 

V • X {(— ) Ta cosnd + Bn slnnml + ( — ) fc cosum + D slnnml) OJ 
n-o " •' 0 ^ n n 

whara V is tha scalar magnatic potantial. Aq and Bn ara 
coofficiants of tha outor potantial. Cn nnd Dq ara 
coafficiants of tho Innar potantial. o Is tha orbital radius. 

N is tha maiifflum ordar of tha modal (usuallg 8). and 9 is tha 
trua anomaly of tho satollito orbit. 


To simplify tha modal, wa assuma tha rafaronca surfaca to ba 
circular, i. a. o ■ constant ■ 6800 hm. To obtain coafficiants A 
through D. wa parform tha appropriata spatial darivativos of V for tha 
magnatic flald compononts. H and Z. and solvo a laast squaros normal 
aquation using tha appropriata MAOSAT vactor magnatomatar obsarvations. 
Tho matrix aquation to ba solvad is: 

X - CM3JL (4) 


whara 

- CZj. Hi* Z2. H2* Zn. 

ara tha maasurad data valuas of tho horixontal and vortical 
fialds H and Z. rospactivaly. and 

Bji 02^ ^2’ * 

ara tha coafficiants for tha potantial. and 


cm 




•<?) 


slum 




-2(f) 

cosm^^ 

-2(f) 

sin©^ 

COE©, sin©, N(— )sln N©, 

i i O 1 

3 

r 

sln©^ 


cosmj^ 


3 





-2(f) 

-(f) 

coam^ 

c 

-2(°-) 

slr.©^ 

*^*****“ * • • « • • 

0 ^ 







Dl 


r 

N(-)cos 

O fll 
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AaiuMinp m>''rN< 'M« ragraition analgtit can ba complatad by tolving 
tha laa«t *«.yuara« rormal aquation; 

2L - Cm'’'‘M3‘1 CMIT ^ (9) 

Tha coafficianta obtainad by thia ragraaaion ara than aubatltutad into 
aquation 3 and tha araociatad darivativaa to calculata fiald 
comp*:nanta. A val ia obtainad for tha acalar potantial and tha fiald 
componanta dua to tha ring currant* 8q daily variation* and any othar 
long wavo4angth fiald which may ba praaant in tha maaauramant 
raaiduala. Thia calculated fiald valvia ia than aubtractad from tha 
obaarvatioea to obtain a fiald valua corv*actad for long wavalangth 
fialda fro^< tha rino currant and tha daily variation. 

It ahoulii bb^ notad that thia procadura ia aimilar to that uaad for 
adjuatmant of cotantial fiald data to a common altituda. By aaauming a 
conatant radiua for o* wa ara in affact adjuating tha componant data to 
that altitu^'i. 

Uhi^a rafvaaaion approach ia apparantly wall~auitad to routina 
corro^^ion of aatallita data, it ia auacaptibla to aliaaing arrora dua 
to r.ia u«a of tha a ^^na fraquanciaa in tha innar and outar functiona. 
Ir.daad.' an aarly formulation of thia mathod for total fiald intanaity 
failitd for thia raaaon. Simulation taata uaing calculatad data ahowad 
that raformulation in tarma of horizontal and vortical fiald componarta 
raliovad thia problam. Thaaa aimulat>..n taata involvad aaauming valiiaa 
for tha coafficianta to aimulata fiald data along a hypothatical orbit. 
Thaaa data wara aubjactad to tha ragraaaion procadura daacribod abova* 
and tha racovarad coafficianta comparod with tha coafficianta uaad to 
aimulata tha data. Tha aatallita orbit radiua for thaaa taata waa 
aaaumad to ba givan by 

R ■ 6800 * lOOainO* (km). 

To aimulata tha fialda dua to tha ring currant and 8q* twalva 
coafficianta wara apacifiad for tha outar functiona and twalva for tha 
innar functiona (N«3). A total of 912 valuaa around tha orbit wara 
calculatad for tha horizontal fiald componant and tha vartical 
componant. Tha raaulta of thia aimulation taat and othara* including 
taata uaing actual MA08AT data ara daacribad in tha naxt aaction. 


BESULia 

Thia aaction praaanta raaulta of aimulation taata of tha ragraaaion 
procadura* followad by a damonatration of tha ragraaaion/harmonic 
filtar on ona orbit of actual MA08AT data. 

Simulationa 

For tha taat uaing aimulatad data* tha racovarad coafficiant valuaa 
ara ahown in Tab la 1 in conpariaon with thoaa aaaumad valuaa uaad in 
ganarating tha aimulatad data. 
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TABLE 1 

Fiald Cotfficivnti (nT*km) 


order 

A_ 

ftllMilltf- 

Viluil 

c 

D 

RttQYtrU 

A B 

Valuos 

(L. 

J) 

1 

>238600 

149900 

-19090 

- 3801 1 

-238600 

149900 

-19890 

- 3801 

2 

191900 

30210 

12680 

-37260 1 

191900 

30210 

12680 

-37260 

3 

-196000 

182800 

19970 

79910 1 

-196000 

182800 

19970 

79910 

4 

126300 

1336 

-17770 

-92660 i 

128300 

1336 

-17770 

-92680 

8 

93490 

189900 

11820 

49430 I 

93490 

16990 

11620 

49430 

6 

-4718 

122400 

-10190 

-4181 1 

-4719 

122400 

-10190 

-4181 

Tho rogrossion process was able to recover the 

assumod 

starting valuos 


of tho coofflcionts q,uitt woll. In foct< tho only orror 1« ono port in 
4710 In tho sixth harmonic. 

To obtain an ostimato of tho dogroo of indopondonco of innor and 
outor torms of tho tamo ordor# this simulation tost was ropoatod* using 
all twolvo harmonics of tho innor functions and tho first harmonic of 
tho outor functions to simulats tho data. Valuos woro rocovorod for 
all twolvo outor coofflcionts. Tho rosults of this tost aro shown in 
Tablo 2. 

TABLE 2 


Fiold Ooofficionts (nT*km) 



A^ 

AimiBlA- 

J 

ViliitA 

C 

D 

RtCflYirtJ 

A B 

_ 0 . 


1 

-238600 

149900 

-19890 

- 3801 \ 

-238600 

149900 

-19890 - 

3801 

2 

191900 

30210 

0 

0 i 

191900 

30210 

0. 1430 -0. 

1904 

3 

-196000 

182800 

0 

0 1 

-156000 

182800 

0.0233 -0. 

1910 

4 

128300 

1336 

0 

0 1 

128300 

1336 

-0.0496 -0. 

0897 

9 

93490 

189900 

0 

0 ! 

93450 

18590 

-0.0164 -0. 

0963 

6 

-4718 

122400 

0 

0 1 

-4719 

122400 

0. 0470 -0. 

0374 


Tho rogrossion procost was ablo to rocovor tho assumod starting valuos 
of tho innor functions without orror« and tho orror of commission in 
tho rocovorod valuos of tho outor functions is loss than 0.2 nT»km. 
For a goocontric radius of 7000 km (highor than tho NA08AT orbit)> tho 
fiold orror corrosponding to this valuo is minusculo. Thus wo can 
conclude tboro is good orthogonality botwoon innor and outor functions 
of tho samo order. This it probably duo to tho uso of a comploto orbit 
(and comploto cycle of oxtornal fiold variation) of data in tho 
rogrossion filter* since those functions differ only by a factor (o/r) 
instead of (r/o>. The eccentricity of the HA08AT orbit provides about 
200 km difforonco in radius botwoon apogee and porigooi so similar 
separation of functions should obtain when processing actual HA08AT 
data. 
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The ••nvitlvltv of tha rogrottion precoot to oliotlng botwoon the 
horoonlco of Pifforont order It llluttrotod tg o tlaiulatlen In whim 
data wat tlwulatod uting outer and inner functient of ordort 1 through 
9. hut cooffic lento were recovered onlg for ordort 1 through 4. Thit 
olMulatlon toot It tumeiarlted In Table 3. 

TABLE 3 

Field Ccofficlentt (nTekit) 


order 

A 

ftityaiJ 
B _ 

Values 

C 

D 

RBCBYtrid 

A B 

i:- - 

_B 

1 

>238400 

149900 

>19890 

> 3801 1 

>244200 

192700 

-13300 

994 

2 

191900 

30210 

12680 

>37260 1 

148300 

31890 

14830 

-39220 

3 

>194000 

182800 

19970 

79910 1 

>198200 

184100 

18710 

77100 

4 

128300 

1336 

>17770 

-92680 1 

126700 

2374 

-19170 

-91900 

9 

93490 

189900 

11820 

49430 i 

91710 

187100 

14100 

49940 

4 

-4718 

122400 

-10190 

-4181 i 

-10890 

119200 

-8919 

-4139 

7 

-93990 

34280 

-3806 

-772 1 





8 

8014 

1019 

898 

-2082 1 





9 

-2178 

2094 

239 

916 i 






At can be teen, aliating duo to truncation can create to«e tignlflcant 
errort In the recovered coeffic lento. Coefficlentt A6 and Cl toew to 
have the largett errort. equivalent to 0. 9 and 0. 6 nT reopectivelg at a 
radiuo of 7000 km. Even though thete errort are negligible, tlnce theg 
are much larger than the errort teen in previout tlmulatlont. theg 
terve to ralte a caution flag. In the regrettlon of actual data, a 
timilar eituation will be pretented, with higher order information 
being pretent in the field meaturementt in the form of the cruttal 
magnetic anemaliet. If thit information aliatet into recovered 
correction coef f ic ientt. then the correction will remove tome of the 
cruttal tignal. Thit thould be tettee further in future ttudiet. 

MAQ8AT Data 

At a demonttration of the regrettion/harmonic filter method, a 
HA08AT orbit from earlg in the mittion wat choten at random The orbit 
choten wat number 1176 on November 9. 1979. The tubtrack of thit orbit 
it thown in Figure S« ttarting over the Siberian arctic, detcending 
through the Indian ocean, and atcending over the eastern Pacific to end 
again in the Ruttian arctic over the Kara Sea. The dittance between 
the end pointt it onlg about 400 km. The field retldualt for the 
horiiontal (X and Y) and the vertical <Z) componentt are plotted in 
Figure 4 with reference to a circular bateline. The field valuot are 
plotted at a tcale of 90 nT per inch radiallg from the circle (potitlve 
outtlde) and one complete revolution of the circle correepondt to one 
orbit period of about 9997 tecondt. Note the pretence of a ehort data 
gap over the eoutheatt Pacific. Thit gap it lets than the ehorteet 
wavelength of the correction functions, however, to no untoward reeulte 
are expected. Thit plot repretentt a slight decimation of the original 
HA08AT data, at onlg everg third point wat telectod for procetting. 
Even at that, there are more than 2000 data pointt reprotented. Alto 
thown in Figure 6 it the satellite altitude, plotted at 100 km par inch 
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ralativa to a circular hatalin# at 400 km. Thara it no obvioua viaual 
corralation hatMaan tha amgnituda oP tha Piald raaiduala and tha 
altituda of tha apacacraPt. 

At with tha timulation ta«ta< tha raParanca radiuc «mb chotan aa a 
conatant at 6800 km. Tha data war# Pirat procaaaad through tha 
ragraaaion analgaia to racovar 4th ordar coaPficianta for both tha 
innar and outar Punctiona. An additional paramatar. Aq. waa addad to 
thta ragraaaion analgaia to ramova conatant biaaaa in tha raaiduala. 
Tha coaPPiciant valuaa ara ahown in Tabla 4. 

TABLE 4 


Fiald CoaPPicianta (nT*km) 


ordar 


Aaaumad 

B. 

■Yaiiiii 

_ J> _ 

PtCflYtriti 

A B 

VilWH 

C 

D 

0 

0 


- 


J -10 

- 



1 

0 

0 

0 

0 

i 21794 

-44320 

-44608 

-31163 

2 

0 

0 

0 

0 

1 7494 

-4276 

-292 

9643 

3 

0 

0 

0 

0 

1 -18949 

2873 

11331 

-13931 

4 

0 

0 

0 

0 

\ 10797 

2428 

9490 

2392 

9 

0 

0 

0 

0 

t -1432 

309 

-3179 

7977 

6 

0 

0 

0 

0 

1 -996 

-1733 

-889 

-2613 

7 

0 

0 

0 

0 

1 -677 

2247 

-376 

1116 

8 

0 

0 

0 

0 

1 -1979 

-1873 

-292 

-661 

Baaed 

on 

thaaa coaPPiciant 

valuaa. 

a corraction Piald waa calculated 


at aach data point and aubtractad Prom tha obaarvad data to giva a 
corractad raaldual. Typical proPllaa oP tha corraction Pialda Por tha 
horitontal and vartical Pialda ara ahown in Figura 7. Nota that tha 
innar and outar corractiona hava aimilar Praguancy contant. aa 
aKpectad. Tha magnituda oP tha corractiona ia typical Por undiaturbad 
ring currant and dally variation Pialda. 

Figuraa 8 and 9 ahow tha horiiontal and vartical Plaid raaiduala. 
raapactivaly t baPora and aPtar corraction. Nota tha complata ramoval 
oP long wavalangth aignala in tha raaiduala. Tha corractad raaiduala 
ara much raducad. but tha abort wavalangth aignala ara praaarvad. 
including tha obvioua aplkaa in tha polar ragiona cauaad by Piald 
alignad curranta. For tha horiiontal componont. tha RMS oP raaiduala 
ia raducad Prom 28.31 nT to 14.47 nT. a raduction oP 24.33 nT RMS. 
whila tha improvamant in tha vartical componant ia 10.99 nT RMS. Prom 
16. 29 to 12. 06 nT RMS. 

cflhfiLuatotffl 

■acauaa oP a lack oP timely data Prom magnatic obaarvator iaa. tha 
praaantly availabla mathod Por medalling tha dally variation Pialda ia 
not practical Per routine correction of MAO BAT vector data. Ragraaaion 
analyaia oP aatallita data, uaing a long wavalangth medal Por tha ring 
currant, waa Pound to yield aa good or batter raaulta aa a more 
complicated modal uaing apharical harmonic ceoPPicianta which wore 
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Punetion* of mognotic indicoo 08 «nd Dot. Thit rogrootion ogproocli lo 
oioro ooionoblo to routino corrocttono bocouoo of tto tndopondonco of 
indicoo dorlvod oo global avoragot from torrcotrial oboorvatorg 
OMgnotograao, Tboroforoi to aoot tho nood for a aoro gonaral 
corroctlen algorltba# and ono wbicb could bo oaplogod on a routino 
baoio to HA08AT vector data by all invootigatoro> a rogrootion/baraenic 
filter procedure wao developed which corrected for both the ring 
current and the dally variation ao well at other long-wavelength 
etternal floldi. 

This procedure ueee the hnown epectrel content of tho etternal field 
current diotributiono to define the order end fora of the 'oatheaetioal 
foraulation« and ueee the obeerved field coaponent reeiduale to 
deteralne the optiaua coefflclente for r.orreetlon of data In that 
paroicjlar HA08AT orbit. In addition to correcting for long-wavelength 
external field varlatlonoi it aleo hao the capability of reducing all 
ebeervatlone to a coaaon altitude on an orbit by orbit baele* without 
prior knowledge of conventional geoaagnetic field indicee. 

Slaulation teete of thie procedure ehow good eeparatlon between 
inner and outer functione of the eaae frequency* and excellent recovery 
of elaulated external field aodel coefflclente. Allaelng due to 
truncation* which will certainly occur in operational uee of the 
procedure* ehould cauee correction errore of leee than 1 nT. Although 
thle le relatively email* it ehould be reevaluated in a more reelletlc 
teet ueing actual data whoee epectrel characterietice ere known. 

Teete of thie procedure on a HA06AT orbit eelected at random ehowe 
vieible improvement In the field component reeiduale* and qualitatively 
good correction for the long wavelength external fielde expected from 
the ring current and the daily variation. Horliontal coaponent 
reeiduale are reduced by 24.33 nT RH8* while the vertical component 
reeiduale are reduced by 10.99 nT RH8. Theee correctione agree with 
the expected aiagnltude of the combined ring current and daily variation 
field. Further evaluation of thie procedure ehould irtclude 
intercomparieone between colinear orbite before and after correction* 
and evaluation of croee-track gradlente between adjacent orbite. 

Thue the objective of thie inveetigation* development of a method 
for routine removal of the meet pereletent external effocte from the 
MA08AT data* hae been eatiefledi deepite c ircumetancee which forced a 
change from the propoeed approach. Not only hae euch a method been 
developed* but it hae been teetod in eimulatlene to examine Ite 
limitatiene and prove feaeiblllty. Ite practical applicability hae 
alee been demenetrated by teete with actual NA08AT data. The harmonic 
regreeeion method alee hae the potential for reduction of the eatellito 
vector data to a common altitude automatically* in addition to 
correcting for long wavelength external fielde. 
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Figure 1. Scalar magnetic residual profiles from 3 colinear POGO passes, at about 
80° W longitude, altitude 470 to 520 km. Note the presence of the 
equatorial electrojet of the daily variation at the magnetic equator, 
10° S latitude (R. Lar ^el, personal communication). 
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Figure 2. Component magnetic fields of the Sq current 

system, December 27, 1964. Note the variation 

with latitude, 0*^ to 40^' North, and the constancy 
with altitude. 
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Scalar magnetic field of the Sq current system on May 24, 1958 
a period of high solar activity. 
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Figure 4. Scalar magnetic residual profiles from 4 collnear P(X»0 passes, at about 82 W longitude, 
altitude 450 to 650 km (R. Langel, personal comiDunlcatlon) . Tbe solid line shows the 
ring current correction calculated by Cain and Davis (1973). Dst indices are from 
Susiura and Poroa 
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Figure 5. Subtrack of MAGSAT orbit 1176, November 5, 1979, for one orbit period. 
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Figure 6. Horizontal and vertical component field residuals 

for MAGSAT orbit 1176, uncorrected for external fields. 
Field magnitude is plotted radially, (positive outward) 
at 50 nT per inch, relative to the circular baseline. 

The orbit starts at the top and proceeds clockwise. Note 
the evidence of large 3rd and Ath order harmonics in the 
residuals. Also shown is the satellite altitude rela- 
tive to a baseline of AOO km, at 100 km per inch. 
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Figure 7. Typical correction profiles for the vertical and 

horizontal field components, outer and inner functions. 
Note the large 3rd and Ath order harmonics and the simi- 
larity in frequency content between the outer and inner 
functions. The peaks near the equator may correspond 
to the equatorial electrojet. The corrections magnitude 
is typical of the undisturbed ring current and the 
dally variation. 
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Figure 8. Horizontal component field residuals, orbit 1176, before 
and after correction. Note the dramatic reduction In 
negative bias In the residuals, and the dlmunltlon of the 
Ath order harmonic behavior. Large peaks near the poles 
are thought to be due to field-aligned currents. The 
RMS of residuals before and after Is 28.31 and 1A.A7 nT 
respectively, for a reduction of 2A.33 RMS. 



Figure 9. Vertical component field residuals, orbit 1176, before 
and after correction. Note the dramatic removal of 3rd 
order harmonic behavior, while preserving the shorter 
wavelength signals. The RMS of residuals before and 
after is 16.29 and 12.06 nT respectively, for an 
improvement of 10.95 nT RMS. 







